BRDTR-331S 


AD-A246  665 

llilfllllillllllllllil! 


TECHNICAL  REPORT  BRL-TR-3315 


BRL 


THE  TENSOR  EQUATION  AX  +  X A  =  G 


MIKE  S  CHE  IDLER 


FEBRUARY  1992 


DTlC 

ELECTS  f 

AR0  2l99i| 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  IS  UNLIMITED. 


U.S.  ARMY  LABORATORY  COMMAND 


BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN  PROVING  GROUND,  MARYLAND 


*  3  « 


■'  '.!  B 


92-04995 

li 


92  2 


NOTICES 


Destroy  this  report  when  it  is  no  longer  needed.  DO  NOT  return  it  to  the  originator. 


Additional  copies  of  this  report  may  be  obtained  from  the  National  Technical  Information  Service, 
U.S.  Department  of  Commerce,  5285  Port  Royal  Road,  Springfield,  VA  22161. 


The  findings  of  this  report  are  not  to  be  construed  as  an  official  Department  of  the  Army  position, 
unless  so  designated  by  other  authorized  documents. 


The  use  of  trade  names  or  manufacturers’  names  in  this  report  does  not  constitute  indorsement 
of  any  commercial  product. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No  0704-0188 


Public  reporting  Durden  for  this  collection  of  information  is  estimated  TO  average  1  *Our  Per  response,  including  the  time  tor  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden.  to  Washington  Headauarters  Services.  Directorate  for  information  Operations  and  Reports.  1215  Jefferson 
Davis  Highway.  Suite  1204,  Arlington,  va  22202-4302.  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0 188),  Washington,  DC  20503 


AGENCY  USE  ONLY  (Leave  blank) 


4.  TITLE  AND  SUBTITLE 


2.  REPORT  DATE 

February  1992 


The  Tensor  Equation  AX  +  XA  *  G 


3.  REPORT  TYPE  AND  OATES  COVERED 

Final  March-Julv  1991 


5.  FUNDING  NUMBERS 

PR:  1L161102AH43 


6.  AUTHOR(S) 


Mike  Scheidler 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADORESS(ES) 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10.  SPONSORING  /  MONITORING 
AGENCY  REPORT  NUMBER 


U.S.  Army  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21005-5066 


BRL-TR-3315 


12«.  DISTRIBUTION /AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  is  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

We  study  the  second-order  tensor  equation  AX  +  XA  =  G  for  symmetric,  positive-definite  A  and  arbitrary  G. 
Motivated  by  applications  in  the  continuum  mechanics  literature,  we  also  examine  some  special  cases  where  G 
depends  on  A  and  another  tensor  H.  For  arbitrary  dimensions,  we  establish  relations  between  the  solutions  X  for 
various  forms  of  G.  These  results,  together  with  Rivlin's  identities  for  tensor  polynomials  in  two  variables,  are 
applied  in  two  and  three  dimensions  to  obtain  explicit  formulas  for  J  in  direct  (component-free)  notation.  The  results 
include  new  formulas  as  well  as  new  derivations  of  previously  known  formulas.  An  application  to  the  kinematics  of 
rigid  motions  is  considered. 


14.  SUBJECT  TERMS 

continuum  mechanics,  kinematics,  strain  rate,  tensors,  linear  algebraic  equations 

17.  SECURITY  CLASSIFICATION 
OF  REPORT 

IS.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

15.  NUMBER  OF  PAGES 

31 


IB.  PRICE  CODE 


NSN  7540-01-280-5500 


Standard  Form  298  {Rev  2-89) 

Prescribed  by  ANSI  Std  239-18 
298  102 


INTENTIONALLY  LEFT  BLANK. 


II 


TABLE  OF  CONTENTS 

Page 


ACKNOWLEDGMENT  .  v 

1.  INTRODUCTION  . 1 

2.  ALGEBRAIC  PRELIMINARIES .  4 

3.  THE  TENSORS  LA  AND  MA  .  6 

4.  SOME  USEFUL  IDENTITIES .  13 

5.  FORMULAS  FOR  THE  THREE-DIMENSIONAL  CASE .  16 

6.  FORMULAS  FOR  THE  TWO-DIMENSIONAL  CASE  .  22 

7.  AN  APPLICATION  TO  THE  KINEMATICS  OF  RIGID  MOTIONS .  25 

8.  REFERENCES .  27 

DISTRIBUTION  LIST .  29 


r 


Intentionally  left  blank. 


iv 


ACKNOWLEDGMENT 


I  would  like  to  thank  Joseph  Santiago  for  his  criticism  of  an  earlier  draft. 


v 


INTENTIONALLY  LEFT  BLANK. 


VI 


1.  INTRODUCTION 


Tensor  equations  of  the  form 


AX  +  XA  =  G 


arise  in  a  variety  of  problems  in  continuum  mechanics;  cf.  Sidoroff  (1978),  Guo  (1984),  Hoger  and 
Carlson  (1984),  Leonov  (1976),  Stickforth  (1982-1983),  Mehrabadi  and  Nemat-Nasser  (1987),  Dienes 
(1979),  Stickforth  and  Wegener  (1988),  Hoger  (1986),  Nemat-Nasser  (1990),  Truesdell  (in  press), 
Wheeler  (1990),  and  Scheidler  (to  be  published)  for  applications  to  kinematics  and  constitutive 
modeling.  Here  the  second-order  tensors  A  and  G  are  regarded  as  known,  and  we  wish  to  solve  (1.1) 
for  the  second-order  tensor  X.  For  the  applications  considered  in  the  aforementioned  references,  A  i? 
symmetric  and  positive-definite.  This  restriction,  which  will  be  imposed  here  also,  guarantees  that  a 
solution  X  of  (1.1)  exists  and  is  unique.  Indeed,  relative  to  any  principal  basis  for  A,  this  solution  is 
given  by  the  simple  component  formula 


a  +  a. 

•  i 


where  ai  are  the  (necessarily  positive)  eigenvalues  of  A.  Of  course,  to  actually  compute  the 
components  of  X  by  means  of  (1.2),  we  must  first  solve  for  the  eigenvalues  and  eigenvectors  of  A. 
This  has  motivated  several  authors  to  seek  explicit  formulas  for  the  tensor  X  directly  in  terms  of  the 
tensors  A  and  G.  An  example  is  the  formula  of  Sidoroff  (1978)  and  Guo  (1984)  for  the  case  where  G 
is  skew  and  the  underlying  vector  space  is  three-dimensional: 


(tf  -  II,)  G  -  (A2G  +  GA2)  . 


Here  I,,,  11^,  and  IIIA  denote  the  principal  invariants  of  A.  Since  equations  of  the  form  (1.3)  are  often 
said  to  be  displayed  in  direct  notation,  we  will  refer  to  formulas  of  this  type  as  direct  formulas  for  X 
or  direct  solutions  of  (1.1). 

Although  the  derivation  of  the  component  formula  (1.2)  is  a  trivial  exercise,  the  derivation  of 
direct  formulas  for  A”  is  a  nontrivial  problem.  Furthermore,  whereas  the  component  formula  (1.2)  is 
valid  for  any  finite-dimensional  vector  space,  the  direct  formulas  become  more  complex  as  the 


dimension  of  the  underlying  vector  space  increases.  For  example,  in  the  two-dimensional  case  the 
solution  of  (1.1)  for  skew  G  is 


X  =  J_  G  ,  (1.4) 

^  A 

which  is  substantially  simpler  than  its  three-dimensional  counterpart  (1.3).  The  formula  (1.4)  is  a 
special  case  of  the  second  of  two  direct  formulas  obtained  by  Hoger  and  Carlson  (1984)  for  arbitrary 
G  in  the  two-dimensional  case.  For  the  three-dimensional  case,  Sidoroff  (1978),  Hoger  and  Carlson 
(1984),  Leonov  (1976)  and  Stickforth  (1982-1983)  have  obtained  direct  solutions  of  (1.1)  for  arbitrary 
G.  For  the  two-dimensional  and  three-dimensional  cases,  Mehrabadi  and  Nemat-Nasser  (1987) 
obtained  direct  formulas  for  X  when  G  =  AH  -  HA  with  H  arbitrary.  For  the  three-dimensional  case 
with  G  (and  thus  X)  skew,  Sidoroff  (1978),  Guo  (1984)  and  Dienes  (1979)  obtained  formulas  for  the 
axial  vector*  of  X  in  terms  of  A  and  the  axial  vector  of  G. 

The  direct  solutions  of  (1.1)  have  been  used  to  obtain  direct  formulas  for  \arious  tensors  of 
interest  in  continuum  mechanics.  Examples  include  formulas  for 

(1)  the  material  time  derivatives  of  the  stretch  and  rotation  tensors  in  terms  of  either  the  velocity 
gradient  (Guo  1984;  Stickforth  and  Wegener  1988)  or  the  stretching  and  spin  tensors 
(Mehrabadi  and  Nemat-Nasser  1987;  Dienes  1979;  Stickforth  and  Wegener  1988;  Hoger 
1986)  or  the  material  time  derivatives  of  the  Cauchy-Green  tensors  (Hoger  and  Carlson 
1984); 

(2)  the  stretching  tensor  in  terms  of  the  Jaumann  rate  of  the  left  Cauchy-Green  tensor  (Sidoroff 
1978;  Leonov  1976); 

(3)  the  derivative  of  the  square  root  of  a  tensor  (Hoger  and  Carlson  1984); 

(4)  the  elastic  and  plastic  spin  tensors  in  terms  of  the  corresponding  stretching  tensors  (Nemat- 
Nasser  1990); 


*The  axial  vector  of  a  skew  tensor  G  is  the  unique  vector  g  such  that  Gv  =  g  x  v  for  every  vector  v,  where  g  x  v  denotes  the 
vector  ctoss  product  of  g  and  v. 
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(5)  the  spin  tensor  in  terms  of  the  rotational  momentum  and  the  Euler  tensor,  for  a  body 
undergoing  a  rigid  motion  (Truesdell,  in  press). 

The  derivatives  of  the  stretch  and  rotation  tensors  with  respect  to  the  deformation  gradient  (Wheeler 
1990;  Scheidler,  to  be  published)  also  satisfy  tensor  equations  of  the  form  (1.1). 

Now  the  tensor  G  on  the  right-hand  side  of  (1.1)  may  depend  on  A  and  some  other  tensor,  say  H, 
which  is  regarded  as  known.  Expressions  for  G  of  the  form 

AH  -  HA,  AHA,  HA  -  AH1,  AH  -  HrA, 

AH  +  HrA  HA  +  AHt,  A2H  +  HA2  (1.5) 

occur  in  the  applications  listed  above;  here  HT  denotes  the  transpose  of  H.  Of  course,  to  obtain  a 
direct  solution  of  (1.1)  in  these  cases  we  can  simply  substitute  the  expression  for  G  in  terms  of  A  and 
H  into  one  of  the  direct  solutions  which  are  valid  for  any  G.  However,  in  many  of  these  cases  there 
exist  direct  formulas  for  X  which  are  simpler  than  (but  necessarily  equivalent  to)  the  formulas  obtained 
by  this  approach.  One  of  the  goals  of  this  paper  is  to  develop  methods  which  yield  these  simpler 
formulas  with  a  minimal  amount  of  computational  effort 

The  paper  is  organized  as  follows.  Section  2  contains  some  preliminary  material.  In  Sections  3 
and  4  we  study  the  general  properties  of  the  solutions  of  the  equation  (1.1)  and  the  equation  obtained 
from  (1. 1)  by  setting  G  =  AH  -  HA,  i.e.,  the  tensor  equation 

AX  +  XA  =AH  -  HA  .  (1.6) 

By  "general"  properties  we  mean  properties  which  are  independent  of  the  dimension  of  the  underlying 
vector  space.  We  show  that  when  G  is  given  in  terms  of  A  and  some  tensor  H  by  one  of  the 
expressions  in  (1.5),  the  solution  of  (1.1)  has  a  simple  expression  in  terms  of  the  solution  of  (1.6). 
Furthermore,  we  show  that  the  solution  of  (1.1)  for  arbitrary  G  can  be  expressed  in  terms  of  the 
solution  of  (1.6)  with  H  =  G. 

In  Sections  5  and  6  we  restrict  attention  to  the  three-dimensional  and  two-dimensional  cases, 
respectively,  since  these  are  the  cases  of  most  interest  in  continuum  mechanics.  We  utilize 
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Rivlin’s  (1955)  identities  for  tensor  polynomials  in  two  variables  to  derive  direct  solutions  of  (1.1)  for 
the  special  case  where  G  is  orthogonal  to  every  polynomial  These  formulas  are  then  used  to 
derive  direct  solutions  of  (1.6)  for  arbitrary  H,  and  those  formulas  are  in  turn  used  to  derive  direct 
solutions  of  (1.1)  for  arbitrary  G.  When  G  is  given  by  one  of  the  expressions  in  (1.5),  simpler  direct 
solutions  of  (1.1)  follow  from  the  identities  in  Section  4  and  our  direct  solutions  of  (1.6).  Several  of 
the  direct  formulas  derived  in  this  paper  have  been  obtained  previously  by  some  of  the  authors 
mentioned  above.  However,  our  derivations  differ  from  theirs,  and  we  derive  some  new  formulas  as 
well.  Our  emphasis  on  the  special  case  (1.6)  of  the  tensor  equation  (1.1)  was  motivated  by  some 
results  in  the  paper  of  Mehrabadi  and  Nemat-Nasser  (1987);  see  the  comments  following  equation 
(4.5)  in  this  paper. 

In  Section  7  we  use  one  of  our  formulas  to  derive  a  new  direct  formula  for  the  spin  tensor  in 
terms  of  the  rotational  momentum  and  the  tensor  of  inertia  for  a  body  undergoing  a  rigid  motion.  In  a 
follow-up  paper  (Scheidler,  to  be  published)  we  use  the  general  properties  and  the  direct  formulas 
obtained  here  to  derive  direct  formulas  for  the  derivatives  of  the  stretch  and  rotation  teasors  with 
respect  to  the  deformation  gradient.  We  apply  those  results  to  the  derivation  of  direct  formulas  for  the 
material  time  derivatives  of  the  stretch  and  rotation  tensors  and  also  to  the  derivation  of  direct 
formulas  for  the  stress  in  a  hyperelastic  material  in  terms  of  the  derivative  of  the  strain  energy  with 
respect  to  the  right  stretch  tensor.  The  latter  problem  was  solved  recently  by  Wheeler  (1990)  for  the 
isotropic  case. 

2.  ALGEBRAIC  PRELIMINARIES 

The  underlying  vector  space  is  a  finite-dimensional  inner  product  space  T'ovcr  the  reals.  By  a 
second-order  tensor  we  mean  a  linear  transformation  from  Vinto  ‘V.  The  vector  space  of  all  second- 
order  tensors  is  denoted  by  Lin.  Sym  and  Slew  are  the  subspaces  of  Lin  consisting  of  symmetric  and 
skew  tensors.  Throughout  this  paper  A  denotes  an  arbitrary  symmetric  positive-  definite  tensor  with 
eigenvalues  a,  >  0,  i  =  1,2, ...,  N,  where  N  is  the  dimension  of  V.  The  identity  tensor  is  denoted  by  /; 
B,  G,  H,  and  X  denote  second-order  tensors  which  are  arbitrary  unless  specified  otherwise.  The  inner 
product  of  B  and  H  is 


B  •  H  :=  tr  ( BJH) . 


(2.1) 
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(2.2) 


sym  H  ■=  1  (H  +  H1)  .  skw  H  :=  1  (H  -  HT)  . 

We  record  the  following  identities  for  later  use: 

2  sym  (AH)  =  AH  +  HT A  =  AH  +  HA  -  2  (skw  H)  A 

=  A  sym  H  +  (sym  H)  A  +  A  skw  H  -  (skw  H)  A  ,  (2.3) 

2  sym  (HA)  =  HA  +  AHJ  =  AH  +  HA  -  2A  skw  H 

=  A  sym  H  +  (sym  H)  A  +  (skw  H)  A  -  A  skw  H  ,  (2.4) 

2  skw  (AH)  =  AH  -  HtA  =  AH  -  HA  +  2  (skw  H)  A 

-  A  sym  H  -  (sym  H)  A  +  A  skw  H  +  (skw  H)  A  ,  (2.5) 

2  skw  (HA)  =  HA  -  AHt  =  HA  -  AH  *  2  A  skw  H 

=  (sym  H)  A  -  A  sym  H  +  A  skw  H  +  (skw  H)  A  .  (2.6) 

A  tensor  H  is  a  polynomial  in  A  if  H  =  cm  Am  +  cm_,  Am~ 1  +  ...  +  c,  A  +  c0  I  for  some 

nonnegativc  integer  m  and  scalars  c0.  c, . cm.  The  subspace  of  Lin  consisting  of  all  polynomials  in 

A  is  denoted  by  •F(A).  The  Caylcy-Hamilton  theorem  implies  that 

T(A)  =  span  (I,  A . A*"') 

=  span  (A,  A2 . A*) 

=  span  (...,  A-2,  A'1,  /,  A,  A2.  ...)  .  (2.7) 

Let  T  (A)1  denote  the  orthogonal  complement  of  T  (A),  i.e„  H  e  9  (A)1  iff  H  is  orthogonal  to  every 
polynomial  in  A.  Then  from  (2.7)  and  (2.1)  it  follows  that 
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(2-8) 


He  P(A)l<=*  tr  (AkH)  =  0,  V  k  e  {0,  1,  ....  N-l) 

<=>  tr  (AkH)  =  0,  Vie  (1,2 . N) 

<=>  tr  (AkH)  =  0,  V  integer  k  . 

In  particular,  if  He?  (A)1  then  H  is  deviatoric,  i.e.,  tr  H  =  0.  Also  note  that  every  skew  tensor 
belongs  to  'P(A)J\ 

By  a  fourth-order  tensor  we  mean  a  linear  transformation  from  Lin  into  Lin.  K  denotes  an 
arbitrary  fourth-order  tensor;  the  image  of  He  Lin  under  K  is  denoted  by  K  [H\.  It  is  easily  verified 
that  the  following  two  conditions  are  equivalent: 

K  [Ht\  =  (K  [H])T\ 

K  [sym  H]  =  sym  K  [H]  and  K  [skw  H]  =  skw  K  [ H J  .  (2.9) 

If  these  conditions  hold  for  each  H  e  Lin  we  say  that  K  is  even.  Similarly,  it  is  easily  verified  that 
the  following  two  conditions  are  equivalent: 

K  [//r]  =  -  (K  [H]f; 

K  (sym  H]  =  skw  K  [H]  and  K  (skw  H]  =  sym  K  [//]  .  (2.10) 

If  these  conditioas  hold  for  every  H  e  Lin,  we  say  that  K  is  odd.  The  transpose  of  K  is  the  unique 
fourth-order  tensor  K  T  with  the  property 

B  •  Kt[H]  =  K  [B]>H  ,  V  A,  He  Lin.  (2.11) 

We  say  that  K  is  symmetric  if  K  T  =  K . 

3.  THE  TENSORS  LA  AND  MA 

In  this  section  we  study  the  properties  of  five  fourth-order  tensors  associated  with  any  symmetric 
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positive-definite  second-order  tensor  A.  The  two  most  important  of  these  are  denoted  by  LA  and  MA 
and  characterized  by  the  conditions 


X  =  La  [G]  <=>  AX  +  XA  =  G  ,  (3.1) 

,Y  =  Ma  [H]  <=>  AX  +  XA  =  AH  -  HA  .  (3.2) 

The  other  fourth-order  tensors  are  denoted  by  NA,  BA  and  CA.  NA  is  introduced  primarily  to  simplify 
the  statement  of  some  of  the  identities  relating  LA  and  MA.  BA  and  CA  are  introduced  to  facilitate  the 
definition  of  the  other  tensors  and  the  derivation  of  their  properties. 

The  fourth-order  tensor  BA  is  defined  by 

Ba  ffl  :=  AX  +  XA  .  (3.3) 

Note  that  BA  is  invertible.  Indeed,  the  condition  BA  [XI  =  G  is  just  the  tensor  equation  (1.1)  which,  as 
pointed  out  in  the  introduction,  has  a  unique  solution  X  for  any  given  G.  We  denote  the  inverse  of  BA 
by  la 

La  :=  (Ba)~'  .  (3.4) 

Then  BA  [X]  =  G  iff  X  =  LA  [C],  which  is  equivalent  to  the  statement  (3.1).  If  we  let  /  denote  the 
fourth-order  identity  tensor  on  Lin,  then* 


which  is  equivalent  to  the  identities 


laba  =  bala  =  i , 


(3.5) 


La  [AH  +  HA)  =  A  La  [H]  +  La  [H]  A  =  H  . 


(3.6) 


*The  product  of  two  fourth-order  tensors  is  defined  to  be  their  composition:  (LA  BA)  [JK]:  =  LA  |B^  [A))  for  each  X  e  Lin. 
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The  fourth-order  tensor  CA  is  defined  by 


CA  [H]  :=  AH  -HA  .  (3.7) 

It  is  easily  verified  that  CA  and  BA  commute;  indeed, 

CA  Ba  =  BaCa  =  CA2 .  (3.8) 

We  define  the  fourth-order  tensor  MA  by 

La  Ca  =  CaLa,  (3.9) 

where  (3.9):  follows  from  prc-multiplying  and  post-multiplying  (3.8),  by  LA  and  using  (3.5).  From 
(3.9)  and  (3.7)  we  see  that 

Ma  [ H ]  =  La\AH  -  HA)  -  A  La  [H]  -  LA  [H\  A  .  (3.10) 

Then  (3.2)  follows  from  (3.1)  (with  G  =  AH  -  HA)  and  (3.10),.  From  (3.9),  (3.8),  and  (3.5)  we 
obtain 

Ba  =  Ba  Ma  =  Ca  =  La  Ca 2  =  Ca2  La  ,  (3.1 1 ) 

which  is  equivalent  to  the  identities 

Ma  \AH  +  HA]=A  Ma  [ff]  +  Ma  [H]  A 
=  AH  -  HA 

=  La  [A2H  -  HA2]  =  La  IH]  -  La  [//]  A2  .  (3.12) 

The  fourth-order  tensor  NA  is  defined  by 

Na  ••=  Ma  Ca  =  Ca  Ma  =  La  (Ca)2  ,  (3.13) 

where  (3.13)2  3  follow  from  (3.9).  From  (3.13)  and  (3.7)  we  see  that 
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NA  \H]  =  Ma  [AH-  HA]  =  A  Ma  [H]  -  Ma  [H]  A 

=  La  [A2H  +  HA2  ~  2  AHA]  .  (3.14) 

By  setting  G  =  A2H  +  HA2  -  2  AHA  in  (3.1)  and  using  (3. 14)3,  it  follows  that 

X  =  Na  IH]  <=>  AX  +  XA  =  A2H  +  HA2  -  2  AHA  .  (3.15) 


The  motivation  for  introducing  the  tensor  NA  is  based  not  so  much  on  the  result  (3.15)  as  on  the 
frequent  occurrence  of  the  expressions  in  (3.14),^. 

It  follows  from  die  results  above  that  the  set  of  fourth-order  tensors 


Sa:={Ba,Ca,La,Ma.Na]  (3.16) 

is  commutative.  Each  tensor  KA  e  SA  also  has  the  following  properties: 

(Ka)T=Ka;  (3.17) 

Ka  [ BH ]  =  B  Ka  [H\  ,  if  AB  =  BA  ;  (3.18) 

Ka  [HB]  =  Ka  [H]  B  ,  if  AB  =  BA  ;  (3.19) 

Kqa(/  [ QHQt ]  =  Q  Ka  [H]  Qt  ,  V  Q  e  Orth  ,  (3.20) 


where  Orth  denotes  the  set  of  orthogonal  tensors  in  Lin.  In  other  words,  KA  [H]  is  an  isotropic 
function  of  A  and  H  which  is  linear  in  H.  Furthermore,  BA,  LA  and  NA  are  even,  whereas  CA  and  MA 
are  odd.  The  easiest  way  to  establish  these  properties  is  to  first  prove  them  for  the  tensors  BA  and  CA, 
and  then  use  the  definitions  (3.4),  (3.9),  and  (3.13)  to  prove  the  properties  for  LA,  MA,  and  NA.  For 
example,  it  is  easily  established  that  (3.18)  holds  for  KA  =  BA  or  CA.  To  prove  (3.18)  for  KA  =  LA, 
use  the  result  BA  [BH]  =  B  BA[H],  apply  LA  to  obtain  BH  =  LA  [B  BA  [W]],  and  then  set  H  =  LA  [A") 
to  obtain  B  LA  [X]  =  LA  [BX],  To  prove  (3.18)  for  KA  =  MA,  use  the  results  for  KA  =  CA  and  LA. 

Ma  [BH]  =  La  [Ca  [BH]]  =  La  [B  Ca  [//]] 

=  B  La  [Ca  [H]]  =  B  Ma[H). 
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The  proof  of  (3.18)  for  KA  =  NA  is  similar. 


In  addition  to  the  expressions  (3.10)  for  MA  in  terms  of  LA,  we  also  have 

Ma  [H]  =  H  -  2  La  [H]  A  =  2  A  La  [H]  -  H 

To  prove  (3.21),,  for  example,  use  (3.10),  (3.6),  and  (3.19)  (with  KA  =  LA  and  B  =  A): 

Ma  [H]  =  La  [AH  -  HA] 

=  La  [AH  +  HA]-  2  La  [HA] 

=  H -2  LA[H]  A 

In  addition  to  the  expressions  (3.14),  2  for  NA  in  terms  of  MA,  we  also  have 

Na  [H]  =  AH  -  HA  -  2  Ma  [H]  A 
=  HA  -  AH  +  2  AMa  [H]. 

To  prove  (3.22),,  for  example,  use  (3.14),,  (3.12),  and  (3.19): 

Na  [H]  =  Ma[AH  -  HA] 

=  Ma  [AH  +  HA]-  2  Ma  [HA] 

=  AH  -  HA  -  2  MA[H]  A  . 

From  (3.21)  we  obtain  the  following  expressions  for  LA  in  terms  of  MA: 

La  [H]  =  1(H  -  Ma  (//])  A-] 

=  i.  A~l  (H  +  Ma  [//]) 

=  1  ( A  'H  +  HA  '  +  A-'  Ma  [ H]  -  Ma  [H]  A-')  . 

4 


(3.21) 


(3.22) 


(3.23) 
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Now  from  (3.14)  it  follows  that 


A~l  Na  [H]  A''  =  MA  [H]  A-1  -  A-1  Ma  [H]  .  (3.24) 

Then  by  (3.23)3  and  (3.24)  we  obtain  the  following  expressions  for  LA  in  terms  of  NA: 

La  [ H ]  =  1  ( A  'H  +  HA -  A'1  Na  [H]  A'1) 

=  1  A”1  {AH  +  HA  -  Na  [tf])  A'1  .  (3.25) 

4 


To  establish  our  next  result  we  need  the  fact  that 

A2H  +  HA2  =  2  AHA  <=>  AH  =  HA  .  (3.26) 

This  can  be  verified  by  taking  components  relative  to  a  principal  basis  for  A: 

A2H  +  HA2  =  2  AHA  <=>  (a?  +  a 2  -  2 a,a;)  H(>  =  0 
<=>  (a,  -  a)2  H,;  =  0 
«  {a,  -  a)  Hy  =  0 
<=>  AH  -  HA  =  0. 

Since  LA  is  nonsingular,  from  (3.7),  (3.10),  (3. 1 2)3,  (3.14),  (3.23),  (3.25),  and  (3.26)  we  conclude  that 
the  following  conditions  are  equivalent: 


AH 
A  La  [H] 
A  Ma  \H] 
Ca[H] 
ma  [H] 
Na\H] 


HA  , 

La  \H]  a  , 

±  MA[H]  A  , 

0. 

0, 

0. 
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La[H]  =  1  A-'H  , 

La[H]  =  1  HA1  . 

i*  [ff]  =  JL  (A~lH  +  HA-1)  .  (3.27) 

4 

Now  suppose  that  H  e  S’ (A).  Then  A//  =  //A,  and  thus  L4  [H]  =  A'1  //.  And  since  A  '  e  !P(A), 
it  follows  that  LA  [H]  e  S’ (A).  Conversely,  if  LA  [ H]  e  S’ (A)  then  by  (3.6)2  it  follows  that  H  e  !P(A). 
Thus  we  have  shown  that 


La  [H]  6  !P  (A)  <=>  H  g  S’ (A)  .  (3.28) 

In  other  words,  if  AX  +  A A  =  W,  then  A"  is  a  polynomial  in  A  iff  W  is  a  polynomial  in  A.  From  (3.6)2 
we  also  see  that 


2  tr  (A*+1  La  [«])  =  tr  (A*  H),  V  integer  A:  .  (3.29) 

Then  from  (3.29)  and  (2.8)  we  conclude  that 

La  [H]  e  T(AY-**He  T (A)x  .  (3.30) 

In  other  words,  if  AX  +  XA  =  H,  then  X  is  orthogonal  to  every  polynomial  in  A  iff  H  is  orthogonal  to 
every  polynomial  in  A.  Since  tr  (A*  (AH  -  HA))  =  0  for  any  integer  k ,  from  (2.8)  with 
H  — »  AH  -  HA  we  conclude  that 


CA  [H]  =  AH  -  HA  e  S’  (A)1  . 
Then  from  (3.10),,  (3.14),,  (3.30),  and  (3.31)  we  see  that 

MA[H],NA[H]e  S’ (A)1. 


(3.31) 


(3.22) 
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In  particular,  for  any  tensor  H  the  solution  X  of  (1.6)  is  orthogonal  to  every  polynomial  in  A.  Of 
course,  since  A  and  H  commute  for  any  H  e  T  (4),  from  (3.27)  we  also  have 

CA  [H]  =  Ma  [H]  =  Na  [H]  =  0,  VHe<P(A).  (3.33) 

We  conclude  this  section  with  a  description  of  our  method  for  generating  direct  formulas  for 
La  [H]  and  MA  [H],  Suppose  that  by  some  means  we  have  obtained  a  direct  formula  for  LA  [G] 
which  is  valid  for  any  6t  f  (A)1  but  not  necessarily  for  other  G.  As  we  will  see  in  Sections  5 
and  6,  such  formulas  are  easily  derived  in  the  two-dimensional  and  three-dimensional  cases.  Since 
AH  -  HA  e  :P  (A)1  for  any  tensor  H,  and  since  MA  [//]  =  LA  [AH  -  HA],  by  setting  G  =  AH  -  HA 
in  our  formula  for  LA  [G]  we  obtain  a  formula  for  MA  [ff]  which  is  valid  for  any  tensor  H.  Then 
from  (3.23)  we  obtain  formulas  for  LA  [H]  which  are  valid  for  any  H.  Alternatively,  we  can  use  the 
relation  NA  [H]  =  MA  [AH  -  HA]  to  obtain  a  formula  for  NA  [//],  and  then  use  (3.25)  to  obtain 
formulas  for  LA  [H]  which  are  valid  for  any  H.  In  these  formulas  for  LA  [H]  we  can,  of  course, 
replace  the  A-1  terms  by  a  polynomial  in  A  via  the  Cayley-Hamilton  theorem. 

4.  SOME  USEFUL  IDENTITIES 

Consider  the  tensor  equation 


AX  +  XA  =  G(A,  H)  ,  (4.1) 

where  G(A,  H)  denotes  one  of  the  expressions  listed  in  (1.5).  In  terms  of  the  notation  introduced  in 
the  previous  section,  the  solution  of  (4.1)  is  X  =  LA  [G(A,  //)];  and  for  the  special  case 
G(A,  H)  =  AH  -  HA  this  solution  can  also  be  written  as  X  -  MA  [H],  As  we  will  see  in  the  next  two 
sections,  the  direct  formulas  for  MA  [H]  (for  arbitrary  H)  are  much  simpler  than  the  direct  formulas  for 
La  [G]  (for  arbitrary  G).  This  should  not  be  too  surprising  in  view  of  the  method  described  in  the 
previous  section  for  generating  these  formulas.  For  the  cases  where  G(A,  H)  is  given  by  one  of  the 
other  expressions  in  (1.5),  the  existence  of  relatively  simple  direct  formulas  for  A'  is  due  to  the  fact 
that  there  are  simple  expressions  for  LA  [G(A,  H)]  in  terms  of  MA  (or  NA),  A  and  H.  The  purpose  of 
this  section  is  to  derive  these  identities. 
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We  begin  with  the  case  G(A,  H)  -  AHA,  that  is,  with  the  tensor  equation 


AX  +XA  =  AHA  , 


(4.2) 


which  is  equivalent  to  the  tensor  equation 


A-'X  +  XA~x  =  H  .  (4.3) 

These  equations  arise  in  the  problem  of  finding  formulas  for  the  material  time  derivative  of  the  right 
stretch  tensor  in  terms  of  the  stretching  tensor  (Mehrabadi  and  Nemat-Nasser  1987;  Hogcr  1986). 

Since  the  solution  of  (4.2)  is  X  =  LA  [AHA],  whereas  the  solution  of  the  equivalent  equation  (4.3)  is 
X  -  La- i  [//],  we  see  that 

La-<  [H]  =  La  [AHA]  =A  La[H]  A  ,  (4.4) 

where  (4.4),  follows  from  (3.18)  and  (3.19)  (with  B  =  A).  By  substituting  the  expressions  for  LA  [H] 
in  (3.23)  and  (3.25)  into  (4.4),  we  obtain  the  identities 

La- .  [H]  =  Ji  A  (H  -  Ma  [//)) 

=  l  (H+Ma  [H])  a 

=  i  (AH  +  HA  -  Na  [//])  .  (4.5) 

4 


In  other  words,  the  solution  X  of  (4.2)  and  (4.3)  is  also  given  by  any  of  the  three  expressions  on  the 
right-hand  side  of  (4.5).  The  second  and  third  of  these  identities  are  equivalent  to  ones  obtained  by 
Mehrabadi  and  Nemat-Nasser  (1987)  under  the  assumption  that  H  is  symmetric.  For  the  three- 
dimensional  case,  they  used  these  identities  and  their  direct  formula  for  JtfA  [H]  to  obtain  direct 
solutions  of  (4.3).* 


*Their  equations  (8.8),  (8. 1 2) — (8. 1 3)  and  (8.16)  correspond  to  equations  (4.3),  (4-5),  and  (4.5),  in  this  paper.  In  their  derivation 
of  (8.13),  and  apparently  of  (8.16)  as  well,  they  utilized  the  symmetry  of  H  (=  2  D  in  their  notation).  However,  these  results, 
as  well  as  their  direct  formulas  for  LA. i  [ff],  are  valid  for  arbitrary  H;  see  also  the  comments  following  equation  (5.26)  in  this 
paper. 
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Compared  with  (4.4),  the  formula  for  MA in  terms  of  MA  is  much  simpler 


Ma-i  =  -  Ma  .  (4.6) 

Indeed,  from  (3.2)  we  see  that  X  =  MA  [H]  iff  A'1  X  +  X  A'1  =  A"1  (-If)  -  ( -H )  A-1  iff 
X  =  Ma-\  [-//],  which  yields  (4.6).  Note  that  the  identities  (4.5)  can  also  be  obtained  by  replacing 
A  with  A-1  in  (3.23)  and  then  using  (4.6)  and  (3.14). 

Now  consider  the  case  G(A,  H)  =  A2H  +  HA 2,  that  is,  the  tensor  equation 

AX  +  XA=  A2H  +  HA2  .  (4.7) 

This  equation  arises  in  the  problem  of  finding  formulas  for  the  Jaumann  rate  of  the  left  stretch  tensor 
in  terms  of  the  stretching  tensor  (Scheidler,  to  be  published).  Alternate  expressions  for  the  solution 
X  -  La  [A2H  +  HA2]  of  (4.7)  follow  from  the  identities 

La  [A2H  +  HA2]  =  A2  La  [ H]  +  La  [H]  A2 

=  2  La  [H]  A2  +  (AH  -  HA) 

=  2  A2  La  [H]  -  (AH  -  HA) 

=  1  (AH  +  HA  +  Na  [H])  .  (4.8) 

2 


(4.8),  follows  from  (3.18)  and  (3.19);  (4.8)2  3  follow  from  (4.8),  and  (3.12);  (4.8)4  follows  from 
(3.14)3,  (4.4),  and  (4.5)3. 

Next,  we  derive  alternate  expressions  for  LA  [G(A,  H)]  when  G(A,  H)  -  HA  -  AHr,  AH  -  HTA, 
AH  +  H1 A  or  HA  +  AHT,  for  applications  of  the  tensor  equation  (4.1)  in  these  cases  see  Guo  (1984) 
and  Scheidler  (to  be  published).  We  claim  that 

La  [HA  -  AHt]  -2  La  [skw  (HA)] 

=  skw  H  -  Ma  [sym  H] 

=  skw  (H  -  Ma  [H]) ,  (4.9) 
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La  [AH  -  Ht A]  =  2  La  [skw  (AH)) 

=  skw  H  +  Ma  [sym  H) 

=  skw  (H  +  Ma  [//])  ,  (4.10) 

La  [AH  +  HtA]  =  2  La  [sym  (AH)] 

=  H  -  2  La  [skw  H]  A 
=  sym  H  +  Ma  [skw  H] 

-  sym  (H  +  Ma  [ff])  ,  (4.11) 

La  (HA  +  AHr\  =  2  La  [sym  (//A)] 

=  H  -  2  A  La  [skw  H] 

=  sym  H  -  Ma  [skw  H] 

-sym  (H-Ma(H)).  (4.12) 

These  identities  follow  easily  from  (2.3M2.6),  (3.6),  (3.10),  (3.18),  (3.19)  and  the  fact  that  MA  is  odd 
(see  (2.10)). 

Wc  conclude  this  section  with  two  identities  which  are  used  in  Scheidler  (to  be  published): 

Ma  [AH  -  HtA]  =  2  Ma  [skw  (4/0) 

=  A  skw  H  -  (skw  H)  A  +  Na  [sym  H]  ,  (4.13) 

Ma  (HA  -  AHt\  =  2  Ma  [skw  (HA)] 

=  A  skw  H  -  (skw  H)  A  -  NA  [sym  H)  .  (4.14) 

These  identities  follow  easily  from  (2.5),  (2.6),  (3.12),  and  (3.14). 

5.  FORMULAS  FOR  THE  THREE-DIMENSIONAL  CASE 

In  this  section  we  derive  direct  formulas  for  the  fourth-order  tensors  LA,  MA,  and  NA  under  the 
assumption  that  the  underlying  vector  space  Vis  three-dimensional.  We  begin  by  listing  some 
identities  for  second-order  tensors  which  will  be  used  in  the  derivation  of  these  formulas.  The 
principal  invariants  of  a  second-order  tensor  B  are  denoted  by  lB,  II£  and  IIIB: 
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(5.1) 


Ig  :=  tr  B,  Iig  :=  1  ((tr  B)2  -  tr  (fl2)],  Illg  :=  det  B  . 

From  the  Cayley-Hanilton  theorem  we  obtain 

A3  =  I*  A2  -  II*  A  +  Ill*  /  ,  (5.2) 

and 

III*  A’1  =  A2  -  I*  A  +  II*  /  .  (5.3) 

Since  the  expression  1  *  /  -  A  occurs  frequently  below,  it  is  convenient  to  introduce  a  special  symbol 
for  this  expression: 

A  :=  I*  l  -  A  .  (5.4) 

The  eigenvalues  of  A  are  a2  +  av  a3  +  av  a,  +  a 2;  in  particular,  A  is  symmetric  pcsitive-dcfinite. 

Also  note  that  (Sidoroff  1978;  Guo  1984;  Stickforth  and  Wegener  1988) 

III*  =  I*  II*  -  III*  =  (a,  +  a2)(a2  +  a,)  (a3  +  a,)  ,  (5.5) 

and  dial  (Guo  1984) 

III*  A"1  =  A2  +  II*  /  .  (5.6) 

The  following  identities  arc  due  to  Rivlin  (1955): 

A2//A2  =  II*  AHA  -  Hi*  ( AH  +  HA) 

+  lA2„  A2  +  a, (A,  H)  A  +  III*  I*,,  I  ,  (5.7) 

AZHA  +  AHA2  =  I*  A //A  -  III*  H 

+  I*„  A2  +  a^A,  H)  A  +  III*  \H  1  ,  (5.8) 

A2//  +  HA2  =  -  AHA  +  I a(AH  +  HA)  -  II*  H 

+  IH  A2  +  a 3(A,  tf)  A  +  a<(A,  tf)  /  ,  (5.9) 
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when' 


a, (A.  H)  =  tr  [( IllA  I  -  UA  A)  H] 

~  11*4  ~  H A  * AH  =  ~  Iw2//  (5.10) 

a 2(A,  H)  =  tr  ((A2  -  I aA)H] 

-  I42//  ~  1 A  U/J  =  ~  1,44//  ’  (5.11) 

o ty(A,  H)  =  -  tr  f(IA  I  -  A)  H] 

=  I4//  -  I4  1  h~  ~  U//  •  (5. 1 2) 

«4(A.  //)  =  tr  [(A2  -  \A  A  +  \\A  D  H) 

~  bvif  ~  I4  I4//  +  H4  1//  =  IH4  I4  '//  (5-13) 


Rivlin's  identities,  as  well  as  (5.2M5.5),,  are  valiu  for  any  second-order  tensors  A  and  H.  However, 
we  will  continue  to  assume  that  A  is  symmetric  posilivc-dc finite. 

Observe  that  if  H  e  'J’(A)1  (sec  (2.8))  then  the  second  line  on  the  right-hand  side  of  Rivlin’s 
identities  (5.7y-(5.9)  vanishes.  In  this  case  we  also  have  the  identity 

A  (AH  +  HA)  A  =  II1A  H,  V  H  e  T(A )x  .  (5.14) 

This  follows  by  substituting  (5.4)  into  the  left-hand  side  of  (5.14),  expanding,  and  then  using  (5.5),, 
(5.8),  and  (5.9).  Similarly,  from  (5.5),,  (5.7),  and  (5.9)  we  obtain  the  identity 

(A2  +  \\A  T)  H  (A2  +  l\A  /)  =  111^  (AH  +  HA),  VHef(A)1.  (5.15) 

In  view  of  (5.6)  we  sec  that  the  identities  (5.14)  and  (5.15)  arc  equivalent. 

Now  suppose  that  G  e  'P(A)±,  in  which  case  LA  (G]  e  TiA)1  also  (see  (3.30)).  If  we  set 
H  =  La  [G]  in  (5.14)  we  obtain 

A  (A  La  [G]  +  La  [G]  A)  A  =  1I1A  La  [G]  . 
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But  since  A  LA  [G]  +  LA  [G)  A  =  G  for  any  tensor  G  (s  (3.6)),  we  have  the  following  simple 
formula  for  LA  (G] : 

III*  La  [G]  =  AgA,  V  G  e  :P(A)X  .  (5.16) 

By  substituting  the  definition  (5.4)  of  A  into  (5.16)  and  expanding,  we  obtain  the  alternate  formula 

III*  La  [G]  =  \\G  -\a  (AG  +  GA)  +  AGA,  V  G  e  T(A)L  .  (5.17) 

By  using  Rivlin's  identity  (5.9)  with  H  =  G,  we  see  that  (5.17)  is  equivalent  to  the  formula 

III*  La  [G|  =  (I*  -  11^)  G  -  (A2G  +  GA1),  V  G  e  'P(A)1  .  (5.18) 

In  view  of  (3.1),  the  formulas  (5.1 6)— (5. 18)  yield  simple  direct  solutions  of  the  tensor  equation  (1.1) 
when  G  e  l’(A)1.  Since  any  skew  tensor  belongs  to  T(A)X,  the  formulas  (5.1 6)— <5. 1 8)  hold  in 
particular  for  any  G  e  Skw.  For  this  special  case  the  formula  (5.18)  was  obtained  by  Sidoroff  (1978) 
and  Guo  (1984).  Their  derivations  employ  the  axial  vector  associated  with  a  skew  tensor  and  thus  do 
not  carry  over  to  th>  more  general  case  where  G  need  only  belong  to  !P(A)X 

Since  MA  [//]  =  LA  \ AH  -  HA\,  and  since  AH  -  HA  e  !P(4)X  for  any  tensor  H  (see  (3.31)),  by 
setting  G  =  AH  -  HA  in  (5.16)  and  (5.17)  we  obtain  the  following  simple  formulas  for  MA  [//]: 

III|  Ma  |//|  ^  A  (AH  -  HA)  A 

=  I‘  (AH  -  HA)  -  lA  (A2H  -  HA2)  +  A2HA  -  AHA2  .  (5.19) 

Furthermore,  since 

AH  -  HA  =  HA  -  Ah  ,  (5.20) 

by  (5.19),  we  also  have 

hi*  ma[H]  =  A  (hA  -  Ah)  A 

=  AhA2-A2HA.  (5.21) 
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The  formula  (5.19)2  was  obtained  by  Mehrabadi  and  Nemat-Nasser  (1987)  by  repeated  applications  of 
the  Cayley-Hamilton  theorem.  We  emphasize  that  (5.19)  and  (5.21)  are  valid  for  any  tensor  H.  In 
view  of  (3.2),  these  formulas  yield  simple  direct  solutions  of  the  tensor  equation  (1.6).  Furthermore, 
in  view  of  (3.1)  and  the  identities  (4.5)u,  (4.9),  (4.10),  (4.11)2-4,  and  (4.12)2^,  the  formulas 
(5.16H5.18),  (5.19)  and  (5.21)  yield  simple  direct  solutions  of  (4.1)  when  G(A,  H)  =  AHA, 

HA  -  AHr,  AH  -  HrA,  AH  +  HrA  or  HA  +  AH1. 

Since  NA  [H\  =  MA  [AH  -  HA],  by  replacing  H  with  AH  -  HA  in  (5.19)]  we  obtain  the  formula 

III*  Na  [H]  =  A  (. A2H  +  HA2  -  2  AHA)  A  .  (5.22) 

By  substituting  (5.4)  into  (5.22),  expanding,  and  then  using  (5.2)  to  reduce  the  A 3  terms,  we  obtain 

IIIA  Na[H]  =  -  2A2HA2  +  21  A  ( A2HA  +  AHA2)  -  2  (I*+  IIA)  AHA 

+  (IA  IIA  +  IIIA)  (AH  +  HA)  -  21a  IIIa  H  .  (5.23) 

Substitution  of  Rivlin’s  identities  (5.7)  and  (5.8)  into  (5.23)  yields 

IIIA  Na  [H]  =  -  4IIa  AHA  +  (IA  IIA  +  3IIIa)  (AH  +  HA)  -  4IA  IIIA  H 

-  2o2(A,  H)  A2  +  205^4 ,  H)  A-  2IIIa  H)  I ,  (5.24) 

where  02  and  CX3  arc  defined  by  (5.11)  and  (5.12),  and 

(^(A,  H)  =  xi  [(IA  A2  +  (IIA  -  IA)  A  -  IIIA  T)  H] 

=  U  ^ a2h  +  -  I4)  I ah  “  mA  \H  .  (5.25) 

By  using  NA  [H]  =  MA  [AH  -  HA],  (5.20)  and  (5.21)  with  H  — »  HA  -  AH,  we  also  obtain 

hi*  Na  [H]  =  A  (A2h  +  hA2  -  2 Ah  A)  A 
=  A3hA  +  AhA3  -  2 A2hA2 

=  -  2  A2hA2  +  1*  (A2hA  +  AhA2)  -2  ii*  AhA  +  in*  (Ah  +  hA)  ,  (5.26) 
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where  the  last  formula  follows  from  the  Cayley-Hamilton  theorem  for  A.  In  view  of  (4.5)3  and  (4.8),, 
the  formulas  (5.22)-(5.26)  yield  direct  solutions  of  the  tensor  equations  (4.2),  (4.3),  and  (4.7).  The 
direct  solution  of  (4.2)  or  (4.3)  obtained  by  substituting  (5.23)  into  (4.5)3  is  equivalent  to  the  first  of 
two  formulas  obtained  by  Mehrabadi  and  Nemat-Nasser  (1987);  cf.  equation  (8.8)  and  the  formula 
preceding  (8.17)  in  their  paper.*  Their  second  formula  is  similar  to  the  one  obtained  by  substituting 
(5.24)  into  (4.5)3;  Rivlin’s  identity  (5.9)  may  be  used  to  obtain  their  result  from  ours  and  vice  versa. 

Now  we  derive  some  formulas  for  LA  [//]  which  arc  valid  for  any  tensor  H.  Substitution  of  (5.23) 
into  the  expression  (3.25)!  for  LA  in  terms  of  NA  yields 

2  III*  La  [H]  =  AHA  -  I a(AH  +  HA)  +  (I*  +  II*)  H 

+  lA  III,,  A']  HA -  in,  (A~'H  +  HA  ' )  .  (5.27) 

This  formula  was  stated  without  proof  by  Leonov  (1976)  and  Stickforth  and  Wegener  (1988).  Leonov 
attributes  (5.27)  to  L.  M.  Zubov;  Stickforth  and  Wegener  refer  the  reader  to  Stickforth  (1982-1983) 
for  a  proof.  Substitution  of  the  formula  (5.3)  for  A'1  into  (5.27)  yields 

2  IIIA  IIIA  La  [H]  =  lA  A2HA 2  -  I*  (A2 HA  +  AHA2) 

+  III*  (A2H  +  HA2)  +  (I*  +  III,)  AHA 

-  ll  II,  (AH  +  HA)  +  (I*  III,,  +  II,  III*)  H  .  (5.28) 

This  formula  is  due  to  Hogcr  and  Carlson  (1984);  they  derived  it  under  the  assumption  that  H  is 
symmetric  and  then  observed  that  the  formula  is  valid  for  any  H.  In  view  of  (3.1),  the  formulas  (5.27) 
and  (5.28)  yield  direct  formulas  for  the  solution  X  of  the  tensor  equation  AX  +  XA  =  H  for  arbitrary 

H 


Observe  that  when  G  e  TiA)1,  the  simple  formulas  (5.16M5.18)  for  LA  [C]  do  not  follow 
immediately  from  (5.27)  or  (5.28).  However,  Hoger  and  Carlson  (1984)  noted  that  substitution  of 
Rivlin’s  identities  (5.7H5.9)  into  (5.28)  yields 


*The  last  term  on  the  right-hand  side  of  their  first  formula  contains  a  misprint;  the  term  should  read  I„  III,, 


b 
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(5.29) 


III*-  La  [H]  =  -  (A2H  +  HA2)  +  (I!  -  II A)  H 

+  p^A,  H)A2  +  p,(A,  H)  A  +  p0(A,  H)  1  ; 

cf.  their  equation  (2.7)  for  expressions  for  the  isotropic  scalar  functions  P,.  They  observed  that  if 
H  e  Skw  then  the  P;  are  zero,  and  consequently,  (5.29)  reduces  to  (5.18)  for  skew  G.  Indeed,  it  is 
clear  from  their  expressions  for  the  Pj(A,  H)  that  they  reduce  to  zero  for  any  He  IP  (A)x,  and  thus 
that  (5.18)  is  a  special  case  of  (5.29).  By  using  a  representation  for  isotropic  tensor  functions  of  two 
symmetric  tensors,  together  with  his  formula  for  Lx  [G]  when  G  e  Skw  (see  (5.18)),  Sidoroff  (1978) 
also  obtained  (5.29).  His  expressions  for  the  P,  are  more  compact  than  those  of  Hoger  and  Carlson; 
however,  Sidoroff’s  expressions  involve  A-'  and  A~2  and  thus  may  be  less  useful  for  some 
applications.  Other  direct  formulas  for  Lt  in  terms  of  A  and  A1  follow  from  (3.23)l>2  and  (5.19)!: 

2  La[H]  =  [H-  -LA  {AH  -  HA)  A]  A"1 
in>i 

=  HA  '  +  __L  A  (H  -  AHA~' )  A 

HI* 

=  A-1  [  H  +  _J_  A  {AH  -HA)  A] 

in* 

=  A~'H+  _L  a  {H  -  A-1  HA)  A  (5.30) 

HI* 

6.  FORMULAS  FOR  THE  TWO-DIMENSIONAL  CASE 

In  this  section  we  derive  direct  formulas  for  the  fourth-order  tensors  LA,  MA,  and  NA  under  the 
assumption  that  the  underlying  vector  space  'F'is  two-dimensional.  We  begin  by  listing  some  identities 
for  second-order  tensors  which  will  be  used  in  the  derivation  of  these  formulas.  The  principal 
invariants  of  a  second-order  tensor  B  are  denoted  by  lB  and  IIB: 

I*:=trfl,  IIS  :=  det  B  .  (6.1) 

From  the  Cayley-Hamilton  theorem  we  obtain 
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II*  A~l  =  I*  /  -  A  . 


(6.2) 


A2=1aA-1IaI, 

The  following  identities  are  due  to  Rivlin  (1955): 

AH  +  HA  =  I*  H  +  l„  A  +  h(A,  H)  I ,  (6.3) 

AHA  =  UaH+IahA-UaIhI  ,  (6.4) 

where 

P,(A,  H)  =  tr  [(A  -  I*  /)  H] 

=  ~  It  hi  =  “  * a~'h  ■  (6.5) 

Observe  that  if  H  e  T(A)L,  then  all  but  the  first  terms  on  the  right-hand  side  of  Rivlin’s  identities 
(6.3)  and  (6.4)  vanish: 


AH  +  HA  =  \aH  ,  AHA  -  \\AH  ,  V  //  e  PiA)1  .  (6.6) 

Now  suppose  that  G  g  V  (A)1-,  in  which  case  LA  [G]  e  <P  (4)x  also.  If  we  set  H  =  L*  [G]  in  (6.6), 
and  use  the  fact  that  ALa  [G]  +  LA  [G]  A  =  G  for  any  tensor  G,  we  obtain  the  following  simple 
formula  for  LA  [G]: 


I aLa[G]  =  G,  VGeffA)1.  (6.7) 

In  particular,  the  formula  (6.7)  holds  for  any  G  e  Skw.  Since  MA  [«]  =  LA  [AH  -  HA],  and  since 
AH  -  HA  e  2>(/l)x  for  any  tensor  H,  by  setting  G  =  AH  -  HA  in  (6.7)  we  obtain  the  following 
simple  formula  for  MA  [H]\ 

\A  Ma  [H]  =  AH  -HA  (6.8) 

This  formula  was  obtained  by  Mehrabadi  and  Nemat-Nasscr  (1987)  using  a  different  method. 

Since  NA  [H]  =  MA  [AH  -  HA],  by  replacing  H  with  AH  -  HA  in  (6.8)  we  obtain  the  formula 

I*  Na  [H]  =  A2H  +  HA 2  -  2 AHA  .  (6.9) 


23 


Substitution  of  (6.2),  into  (6.9)  yields 


lANA[H]  =  -2  AHA  +  1A  (AH  +  HA)  -  2  II AH  .  (6.10) 

Recall  that  X  =  LA- 1  [H]  is  the  solution  of  the  tensor  equations  (4.2)  and  (4.3).  As  an  application  of 
the  fonnula  (6.10)  for  NA  [H],  observe  that  substitution  of  (6.10)  into  the  identity  (4.5)3  yields  the 
formulas 

2  lA  La->  [tf]  =  AHA  +  IIA  H 

-  2  +  lyw  A  -  UA1H  I .  (6.11) 

The  second  fonnula  follows  from  the  first  by  Rivlin’s  identity  (6.4).  A  formula  equivalent  to  (6.11), 
was  obtained  by  Hoger  (1986)  using  a  different  method;  cf.  equation  (4.5)  and  the  equation  preceding 
(4.3)  in  her  paper. 

Finally,  we  derive  some  formulas  for  LA  [H]  which  are  valid  for  any  tensor  H.  Substitution  of 
(6.10)  into  (3.25)2  yields 

2  I*  La  [H]  =  H  +  HA  A-,/fA_I  .  (6.12) 

This  also  follows  from  (6.11),  with  A  -*  A'1  and  the  identities  I*-i  =  IA/1IA  and  II*-i  =  1/TI*.  Recall 
that  A-1  can  be  expressed  in  terms  of  A  by  the  simple  formula  (6.2)2.  By  substituting  (6.2)2  into 
(6.12)  and  expanding,  we  obtain 

2  Ix  II*  La  [H]  =  AHA  -  I*  (AH  +  HA)  +  (I*  +  II*)  H  .  (6.13) 

Substitution  of  Rivlin’s  identities  (6.3)  and  (6.4)  into  (6.13)  yields 

2  lA  IIA  La[H]  =  2  UaH  +  P,(A,  H)  A  +  ^(A,  H)  I .  (6.14) 

where  P,  is  given  by  (6.5),  and 

P0(A,  H)  =  tr  [  ((tf  -  nA)  /  -  IA  A)  H] 

=  (£-«*)  Itf-lxW-  (615) 
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The  formulas  (6.13)  and  (6.14)  were  obtained  by  Hoger  and  Carlson  (1984). 


7.  AN  APPLICATION  TO  THE  KINEMATICS  OF  RIGID  MOTIONS 

The  material  in  this  section,  with  the  exception  of  the  last  equation,  is  taken  from  §  1.10  of  the 
textbook  by  Truesdell  (in  press).  Let  W  denote  the  spin  tensor.  Let  M,  E  and  J  denote  the  rotational 
momentum,  the  Euler  tensor  and  the  tensor  of  inertia  of  a  body  relative  to  the  (possibly  time  - 
dependent)  point  x0.  J  and  E  are  symmetric  positive-definite  tensors  which  are  related  as  follows: 

7  =  (tr  E)l-E.  (7.1) 

If  the  body  is  undergoing  a  rigid  motion,  and  if  x0  is  either  the  location  of  the  center  of  mass  of  the 
body  or  a  fixed  point  of  the  body,  then  the  rotational  momentum  relative  to  x0  is  related  to  the  spin 
tensor  and  the  Euler  tensor  relative  to  x0  by  the  formula 

M  =  -  ( EW  +  WE)  .  (7.2) 

Truesdell  (in  press)  utilized  Guo’s  formula  (see  (5.18))  to  solve  (7.2)  for  the  spin  tensor  in  terms  of 
the  rotational  momentum  and  the  Euler  tensor: 


(7.3) 


The  purpose  of  this  section  is  to  point  out  that  W  has  a  simpler  expression  in  terms  of  the  tensor  of 
inertia.  Indeed,  since  the  solution  of  (7.2)  is  W  =  LE  [-M],  by  setting  A  =  E  and  G  =  -M  in  (5.4), 
(5.5)  and  (5.16)  and  using  (7.1),  we  find  that 

W  *  -  _L  JMJ  .  (7.4) 

III, 
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